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Abstract 

We compute the one-loop unitarity corrections 0(p e ) from K + — > 7r + 7r + 7r~ to 
K + — > 7r + 77 and we find that they are relevant, increasing the leading order pre- 
diction for the width in a 30 — 40%. The contributions of local 0(p e ) amplitudes, 
generated by vector meson exchange, are discussed in several models and we con- 
clude that the vector resonance contribution should be negligible compared to the 
unitarity corrections. 
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1 Introduction 



The phenomenology of radiative non-leptonic kaon decays provides crucial tests for 
the ability of Chiral Perturbation Theory (xPT) |IJ to explain weak low-energy processes 
and interesting possibilities to study CP violation in these channels. xPT is a natural 
framework that embodies together an effective theory (satisfying the basic chiral symmetry 
of QCD) and a perturbative Feynman-Dyson expansion. Its success in the study of 
radiative non-leptonic kaon decays has been remarkable (see 0, |3|] and references therein). 
Of course there are still open problems, but upcoming experiments should improve our 
phenomenological knowledge of these decays, in particular on K — > ^77, K — > tt£ + £~, 
K L ^£+£- or K L ^i£ + £~. 

Kl — > 7r°77 is a very interesting channel by itself as a xPT test and also in order to 
establish the relative role of the CP conserving contribution to Kl — > ir°e + e~ versus the 
CP violating contributions 0, |3|, [4], |f| ^| [7], || . K L — > 7r°77 has no 0{p 2 ) tree level con- 
tribution since the external particles are neutral. For this same reason there are no 0(p 4 ) 
counterterms and the chiral meson loops are finite at this order ||. Consequently xPT 
gives an unambiguous prediction. Furthermore only an helicity suppressed amplitude for 
Kl — * ir°e + e~ is generated, which is small compared to the CP violating ones. At higher 
orders in xPT a new invariant amplitude, generating an unsuppressed helicity amplitude 
to Kl — > 7T°e + e~, appears. As we will see in detail later there is a specific kinematical 
region in the diphoton invariant mass spectra of the K — > rcy-/ processes where only this 
last amplitude contributes. This is the region where both photons are nearly collinear, 
i.e. at small diphoton invariant mass. The experimental results showed |10| that while 
the 0(p A ) predicted spectrum looked in good agreement with the experiment, finding no 
contribution at small diphoton invariant mass, the predicted branching ratio was under- 
estimated by a factor slightly bigger than two, two sigmas away from the experimental 
result. 

This fact prompted to several authors to consider higher chiral order corrections: i) 
Vector Meson Exchange contributions to the local 0{p & ) amplitude (7], ||, ii) Unitarity 
corrections from Kl — > 7r°7r + 7r~ [O, |T^] and iii) Complete unitarization of the tctc in- 
termediate states through a Khuri-Treiman treatment and inclusion of the experimental 
77 — > 7T°7r amplitude [[13]. It has to be emphasized that these higher order corrections 
do not represent a complete chiral order contribution. The size of the first correction is 
controversial 0, § , however a large contribution by itself is excluded by the experimental 
spectrum. The second contribution shows that unitarity corrections of Kl — > 7r°7r + 7r~ to 
Kl — > 7r°77 increase 20 — 30% the leading order amplitude and then have to be taken 
into account [[11], [T^| . When both contributions are added it was possible to fit both the 



width and the spectrum [|12|]. The third correction increases the branching ratio by an 
additional 10 - 20%. 
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The charged channel K + — > 7r + 77 is experimentally less known (only an upper limit 
on the branching ratio exists [0), but measurements with good precision are foreseen 
in the near future ]3], [0|. The leading one-loop 0{p i ) result has been computed || 
and depends upon unknown weak local amplitudes which however could give important 
informations on vector meson weak interactions. Indeed while at 0(p 4 ) it is well known 
that the finite part of the counterterms in the strong xPT lagrangian is saturated by 
the spectrum of lightest resonances (vector, axial-vector, scalar and pseudoscalar) fl6|l , 
the situation in the 0(p 4 ) weak lagrangian is less clear due to our ignorance on the weak 
couplings of vector mesons, however there are predictive models to test [ |17jj . 

We have evaluated the K + — > 7r + 7r + 7r~ unitarity contribution to K + — > 7r + 77 with the 
same procedure used to compute the K L — > 7r°7r + 7r~ unitarity corrections to K L — > 7r°77 



TTjl . The dispersive contribution is unambiguously computed up to a polynomial piece 
which is reabsorbed (together with the loop divergences) in the counterterms. Moreover 
various vector dominance exchange models are studied to saturate the bulk of the coun- 
terterms as in Kl — > 7r°77 . Then we can give a definite prediction for the kinematical 
region of nearly collinear photons where the unknown contribution of the 0(p 4 ) local am- 
plitude is negligible. The experimental study of K + — > 7r + 77 in the DA$NE ^-factory 
and other experiments [[0| is going to improve the phenomenology of K — » Tcy-f and 
will allow us to analyze both channels in a correlated way. 



2 K — > 7T77 amplitudes 

The general amplitude for K — > Tc-f-f is given by 

M ( K(p) -> n(p 3 ) j( qi , e x ) j(q 2 , t 2 ) ) = e lp e 2u M^(p, q u q 2 ) (1) 
where ei,e 2 are the photon polarizations, and M^ v has four invariant amplitudes 

= ((£ q\ - Ql ■ q 29 n + 2 -^A {-p.q ip .q 2 <f> - q x . q 2 p»tf 

mis m, 



+ p- qi q 2 p v + p ■ q2P l "q\ 



+ e^q lp q 2a + [ e^»° (p -q 2 q lp + p- q X q 2p )p„ 

m K F m K F H 



(2) 



+ (p* e ™fh + p v e^)p a q lf3 q 2 



where 



y = , z = fa +f ) 2 . ( 3 ) 

m K m K 

The physical region in the adimensional variables y and z is given by : 

< \y\ < l -\ l/2 (l,rl,z) , < ,2< (l-rv) 2 , (4) 
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with 



X(a,b,c) = a 2 + b 2 + c 2 - 2 (ab + ac + be) 



(5) 



m K 



Note that the invariant amplitudes A(z,y), B(z,y) and C(z,y) have to be symmetric 
under the interchange of q\ and q2 as required by Bose symmetry, while D(z,y) is anti- 
symmetric. In the limit where CP is conserved the amplitudes A and B contribute only 
to Kl — > 7r°77 , while C and D only contribute to Kg — > 7r°77. In K + — > 7r + 77 all of 
them are involved. 

Using the definitions (@J^) the double differential rate for unpolarized photons is given 

by 

d 2 T 



dydz 



m K 
2V 



(| A + B\ 2 + \C\ 



(6) 



y 



B\ 



D 



The processes K — > 7177 have no tree level 0(p 2 ) contribution. At 0(p 4 ) the amplitudes 
-B(z, y) and ^(z, y) are still zero since there are not enough powers of momenta to generate 
the gauge structure, and therefore their leading contribution is 0{p & ). As can be seen 
from (|6]) only the B and D terms contribute for small z (the invariant amplitudes are 
regular in the small y, z region). The antisymmetric character of the D(z,y) amplitude 
under the interchange of q\ and q2 means effectively that while its leading contribution is 
0(p 6 ) this only can come from a finite loop calculation because the leading counterterms 
for the D amplitude are 0(p 8 ). However also this loop contribution is helicity suppressed 
compared to the B term. As shown in a similar situation in the electric Direct Emission 
of Kl — > 7T + 7r~7 [ 18 1 , this antisymmetric 0(p 6 ) loop contribution might be smaller than 
the local 0(p 8 ) contribution. 



3 K + -> 7T+77 at 0{p l 



The leading AI = 1/2 0(p 4 ) A(z,y) and C(z,y) amplitudes for K~ 
already been computed in 0. We review them here. 
The A amplitude reads 

G 8 m 2 K a en 



7T+77 have 



2irz 



{z+l-ri)F{-) + (z + r 2 -l)F(z) 



cz 



(7) 



Here G$ is the effective weak coupling constant determined from K — > nn decays at 0(p 2 



9.2 x 10- 6 GeV~ 2 



(8) 
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and the Fix] function is defined as 



F(x) 



^ arcsm [■\~) 



(3{x) 



1 + 



'i-i 

a; 



In 2 



i + 



7T 



2m In 



l-P(x) 
1 + 



x < A 



x > 4 



(9) 



In © the pion loop contribution F(z/r%) dominates by far over the kaon loop amplitude 
with F(z). The loop results are finite. However as we have already commented xPT 
allows an 0(p 4 ) scale independent local contribution that in (^) is parameterized by 

. 128vr 2 



-[3(L 9 + L 10 )+N U -N 15 -2N : 



18 



(10) 



that is a quantity 0(1). The L 9 and L 10 are the local 0(p A ) strong couplings and N u , jV 15 
and iV 18 are 0(p A ) weak couplings, still not fixed by the phenomenology, and that can be 
only computed in a model dependent way [ 17H . The Weak Deformation Model (WDM) 
H predicts c = 0, while naive factorization in the Factorization Model (FM) |19|, [HJ gives 
c = —2.3. In these models, due to the cancellation in the vector meson contribution in 
(|i0|), the role of axial mesons could be relevant [ |17|] . 
The 0(p 4 ) contribution to the C(z,y) amplitude is 

2 + rl 



C(z) 



G 8 m K a e 

71 



z - K 



r n + ir n - 



m K 



where r v = m^/rriK and T n ° = T(n° — > 77). This amplitude is generated by the Wess- 
Zumino-Witten functional [OT (7r°, 77) — > 77 through the sequence K + — > 7i + (7i°,i]) — > 
7r + 77. This contribution amounts roughly to less than 10% in the total width. 



4 0(p 6 ) local amplitudes for K + 



7T + 77 



At 0(p 6 ) there are only four independent Lorentz invariant local amplitudes contribut- 
ing to K + — > 7r + 77 : 



F^F^d'K+dxir- 



m 2 K F^F^K + n 



F, v F^d x K + d\- 



d a F^d a F^K + Ti- 



(12) 



where the last one has no analogous in Kl — > 7r°77. In general the couplings of these 
operators are not directly related with the ones in Kl — > 7r°77 due to the fact that the 
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electric charge matrix does not conmute with the generators of the electrically charged 
field. As we shall see, in specific models, relations between the two channels can be found. 
All of these operators contribute to A(z,y) but only the first one in (|T2|) gives a B(z,y) 
amplitude. Loop divergences at 0(p 6 ) are absorbed in the counterterm coefficients, that 
thus renormalized are finite. Chiral dimensional analysis [|IJ tells us that their contri- 
butions are suppressed compared to 0(p 4 ) by a factor m 2 K / (Ati F^) 2 ~ 0.2. Nevertheless 
vector meson exchange was found to enhance this up to m 2 K /m\r ~ 0.4 [16|]. Thus we 
try to estimate the contributions of the lightest resonances, i.e. vector mesons, and we 
assume that heavier resonances and non-resonant contributions give smaller corrections. 
This picture seems well verified in the strong coupling constants at 0(p A ) [16], and it is 
likely to apply to the weak couplings. 

We are going to consider here the contribution of vector meson resonances to the weak 
0(p 6 ) counterterms for K + — > 7T + 77 (a more complete discussion and an extension to 
Kl — > 7r°77 is given in [^TJ) and we find that only the terms in ( |12"D with derivatives on 
the meson fields can be generated by vector meson exchange. 

Following the study of K L — > 7r°77 in [|J we assume that the contribution to the local 
amplitudes for K + — > 7r + 77 is dominated by vector meson resonances and define their 
contribution with an adimensional parameter ay generated by the first term in (|T2|) as 

a v = ~ — 2 lim n B v(z) (13) 

2G 8 m K a em *-o 

where By(z) is the vector resonance contribution to the B amplitude. The A amplitude 
generated by vector exchange {Ay) gets contributions from the first and second structure 
in (12). However as seen in || if we assume that these local amplitudes are generated 
through strong resonance exchange supplemented with a weak transition in the external 
legs, those two contributions are related and can be written in terms of the ay parameter 
defined in (0) as 

A v = a y (3 + r w - z) (14) 

7T 



In H two different vector contributions to ay in (|13T ) have been proposed : the first one 
amounts for the weak counterterms generated by a strong vector resonance exchange with 
a weak transition in an external leg (aff*), the second is generated by vector resonance 
exchange between a direct weak vector-pseudoscalar-photon (VP7) vertex and a strong 
one (ay r ). While the first can be computed in a model independent way, the generation 
of direct VP7 weak vertices is still poorly known and therefore only models can be used. 
The external weak transition for K + — > 7r + 77 gives 

= _ 128^,4. = _ 08 

where m v = m p and the strong V^P7 coupling \hy \ — (3.7 ± 0.3) x 10~ 2 is defined by 

£(KP 7 ) = h v e, vpu ( V { u u , /f }) , (16) 
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with 




U = uu , D^U = dJJ - ir^U + im,, , (17) 



U = exp 

and m(0) is an element of the coset space SU(3)l <%> SU(3)r/ SU(3)y parameterized in 
terms of the Goldstone fields fa i — 1, ...8, F^ L are the strength field tensors associated to 
the external r M and £ M fields, is the nonet of vector meson resonances, F n ~ 93 MeV 
is the decay constant of pion and, finally, (A) = Tr(A). In our case with two photon 
external fields = r M = eQA^ [] . 

The model dependent contribution to ay from direct weak vertices in the WDM and 
FM quoted above is : 

dir n ex t 

a V,WDM — a V J 

(18) 

dir _ _OL, „ext 

U V,FM ~ ^tvpUy , 

where in the FM hp is the unknown fudge factor that is not fixed by the model and 
satisfies < kp < 1. Naive factorization predicts kp = 1 . 
By adding ( |15|) and fll8f) we see that ay = a e y l + a* r is 



a V , WDM — U , 

(19) 

«v,fm = {l — 2kp)ay t 

We note that for the allowed range of values of kp in the FM both models agree in 
predicting a very small vector meson contribution to 0(p 6 ) in K + — > 7r + 77 and therefore 
the important conclusion of our exercise is that the vector contribution to the local B(z, y) 
amplitude can be neglected. If, as it is reasonable to assume, other resonance interchange 
corrections are even smaller, we can conclude that the small z region of K + — > 7r + 77 is 
completely predictable at 0(p 6 ) through chiral loops and it is likely to be dominated by 
the unitarity corrections of K + — ► 7r + 7r + 7r~ to K + — > 7r + 77 . Indeed in that region the 
0(p 6 ) is dominant and our ignorance on the c amplitude is irrelevant. 

1 For a full discussion of the definitions and notations see |EgL II. 



6 



5 Unitarity corrections of K + —> 7T + 7r + n to K + — > 

7T + 77 

The amplitude for the process K(p) — > 7r(pi)7r(p2)^(p3) can be expanded in powers of 
the Dalitz plot variables 

x = y = 

7T 7T 

where S{ = (p — j9«) 2 for i = 1, 2, 3 , s = (si + S2 + S3)/3 and the subscript 3 indicates 
the odd pion. For the decay K + {p) — > Tr+Qo^Tr+Qt^^fe) the isospin decomposition, 
neglecting the small phase shifts and up to quadratic terms, is written as |22], ^3| 

_> tt+tt+tt-) = 2 a x - a 3 + (A - i/? 3 + v^s)^ (21) 

-2(Ci + Ca)(^ 2 + ^) " (6 + & - " ^) , 

where the subscripts 1 and 3 refer to A J = 1/2,3/2 transitions respectively, and the 
coefficients in (]2l|) have been fitted to the data [p3~| . The 0(p 2 ) amplitude 



A&>(K+ -> tt+tt+tt-) = G 8 m^ g - r 2 y ) (22) 

with the value of Gg in (H) underestimates by 20-30% the experimental linear slopes for 
K + — > tt+tt+tt - . At the next chiral order, the experimental linear and quadratic slopes 
in (|2lD are recovered (with predictive power too) f23fl . Since the 0(p 4 ) loop contribution 



(0) to K + — > tt+77 , is generated by (|22|) , it seems natural to try to include all the 
contributions to K + — > Tr+77 generated by the experimental slopes in Q2"ID, i.e. we 
evaluate the 0(p 6 ) contributions to K + — > Tr+77 induced by the 0(p 4 ) corrections to 
K + — > tt+tt+tt - . This can be done similarly to the case of Kl — ► tt°77 in |TT| , [T2 



(PH) is considered as an effective i^+ — > tt+tt+tt - chiral vertex, the kinematical variables 
are replaced by the appropriate covariant derivatives, the QED scalar vertices are added 
through minimal coupling and then the usual Feynman diagrams approach can be used. 
There are 13 Feynman diagrams that have two different topologies. There is a subset of 9 
diagrams where one or both photons are radiated by the external legs and therefore do not 
give and absorptive contribution. The sum of these bremsstrahlung-like diagrams is not 
gauge invariant and cancels with an analogous term from the 4 remaining diagrams (two of 
which have also an absorptive part) and generate the A(z, y) and B(z, y) amplitudes. This 
loop result for the A(z, y) and B(z, y) amplitudes is divergent and needs to be regularized 
and renormalized. We have used dimensional regularization and divergences are absorbed 
by the four counterterms discussed earlier. 

Alternatively one could use a subtracted dispersion relation. While the absorptive 
contribution is uniquely determined, the dispersive one, due to the presence of subtraction 
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constants, can be computed only up to a polynomial. This is precisely the same ambiguity 
that happens in the effective vertex method that we have used where the finite part of 
the counterterms is a priori unknown. 

We assume that the dominant contribution should come from the non-polynomial am- 
plitude generated by the cut of the two-pion intermediate state. Otherwise contributions 
generated by a vanishing on-shell K — > 3tt amplitude can be reabsorbed in the unknown 
polynomial amplitude. 

The final result for the A and B amplitudes for K + — > 7r + 77 in the MS subtraction 
scheme are 



^7 



2(2 ttl - a 3 ) + (l + ^2 - ^) (A - \fc + vfyj 



3r 



4 (4Ci 



7T 7T/ 

2 ' 



Hi 



^)+9«-^)-F - 



12 



1 + 6rf) In 



TO, 



+ 



-1(9^-5-3(1 + 3^-,)) |f 



z \12 

Az 



z 



12 



+ R 



z 



1 / 2r; 
+ 2 1 + V 



2 



+ (1 



* + 24 



72r^ 



—2 



36 



2rj 



24z 



(z + r 2 n - l)-F(z) - c + 2r^i + 2r] 2 

z 



(23) 



8 



B(z,y) 



a* 



7T 3r, 



(4<i + (i, 




z \r 



+ G 8 m 2 K r] 3 } , 
where the -F(x) function has been defined in (|9|) and 



(24) 



R(x) 



1,2 2 [A 7 • /v 7 ^ 
£ + — 1 arcsm(-^-) 



6" 



1 
5 



x V x 



2 | /9(g) 
a; ~ r a; 



In 



1 - 
1 + 



+ Z7T 



x < 4 



x > 4 



(25) 



with /?(x) defined in (g) (F(z) ~ - ^ and R{z) 



m ioiz 







In 



the rji 



1, 2, 3 stand for the unknown polynomial contribution we have 



spoken before. Though, in principle, four different subtraction constants could appear, 
these are the only ones that are necessary in order to absorb the loop divergences generated 
by this unitarity correction. We emphasize that these local amplitudes are not generated 
by the vector mesons (already taken into account by ay) and therefore are expected to 
be suppressed by m 2 K /K^, with A x ~ AnF^ over the previous chiral order. If as a naive 
chiral dimensional analysis we choose as coefficients of the suppression factor the factor 
accompanying In (m 2 /fi 2 ) (for definiteness) we get 



12- 



m 



K 



A 



2.06 



7 m 



K 



5 A 



0.24 



V3 = ~ 



3 m 



K 



2 A 



-0.26 



(26) 



We note the big numerical value for rji, but looking into the A amplitude fl23|) we see that 
r)i is suppressed by m^/m 2 ^. The order of magnitude we get indeed is the one expected by 

3/2 coefficients of the 0(p 2 ) 



2jp we only have included the AI 



chiral counting. In 

contribution of K + — > 7t + 7t + tt~ because the bigger errors in the determination of these 
coefficients in the 0(p 4 ) amplitudes. 

We notice that at this order while the B amplitude is y-independent the A amplitude 
gets its first y dependence. 

As we have seen the vector meson contribution to the B amplitude is likely to be very 
small. If we use the same definition (JT3) for our B amplitude (El) we get 



4Cl + 6 

,2 4 



hi 



18G 8 mK V ^ 



;V3 



(27) 



The first term amounts to ~ 0.52 if fi = m p is taken. Then from ( E6|) we can conclude 
that due to the expected vanishing of the vector meson contributions, the polynomial 
non-resonant amplitude, though small compared with the unitarity corrections, could be 
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in principle tested in this channel at low z. We note that this situation is different to the 
Kl — > 7r°77 case where the vector meson contribution seems to be more relevant. 

We have not included the AJ = 3/2 couplings in the 0(p 6 ) contribution because of the 
big errors in the experimental fit. But we should stress that inputting those amplitudes 
could modify the low-z behaviour of the B amplitude even a 30%. 

The argument about the relevance of the correction to Ki — > 7r°77 coming from 
the inclusion of the experimental 77 — > 7r°7r° amplitude f|U|] is weakened here due to the 
relative suppression of K + — > 7r + 7r°7r° compared to K + — > 7r + 7r + 7r _ . 



6 K + —> 7r + 77 : branching ratio and spectra 

We can now proceed to study the numerical results for K + — ► 7r + 77 taking into 
account the 0(p A ) and our 0(p 6 ) loop evaluation of the unitarity corrections. As we have 
shown the 0(p 6 ) vector resonance dominated local contributions are negligible. In our 
numerical discussion we will input rji — 0, i — 1, 2, 3 and u = m p . 

At present there is only an upper bound for By(K + — > 71+77) which depends on the 
shape of the spectrum P^j 

Br(K + -> tt+77 ) | e:E p < 1.5 x 10" 4 (xPT amplitude) , 

(28) 

Br(K + — ► 7r + 77 ) | e:E p < 1.0 x 10 -6 (constant amplitude) 

The uncertainty in the theoretical prediction is dominated by the unknown OQo 4 ) coun- 
terterm generated amplitude c in (0). In Fig. 1 we show Br(K + — > 7r + 77) as a function 
of c with and without the 0(p 6 ) corrections that we have computed. We remind that 
WDM predicts c = while naive FM gives c = —2.3 with the following results : 

Br{K+ -> 7T+77 ) \ WDM = 7.24 x lO" 7 , 

(29) 

Br{K + -»■ 7T+77 ) | nFM = 6.20 x 10" 7 . 

When comparing with the 0(p 4 ) predictions || we find that the unitarity corrections 
increase around 30 — 40% the branching ratio. In Fig. 2 we show the z-distribution at 
0(p 4 ) and with our 0(p 6 ) unitarity correction for c = 0. As can be seen the correction is 
noticeable. 

The uncertainty on the c amplitude also translates into the spectra. In Fig. 3 we show 
the spectrum of the invariant mass of the two photons for the two values of c predicted by 
the WDM and naive FM. We notice that the main dependence on c arises in the z-region 
where there is the bulk of the absorptive contribution. 
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Finally in Fig. 4 we show the spectrum in the "asymmetric" y variable for c = 0. At 
this order this spectrum has much less structure than the one in the di-photon invariant 
mass. 

7 Conclusion 

The situation of the chiral prediction for Kl — > 7r°77 and the expected experimental 
measurement of Y(K + — > 7r + 77) at DA$NE and BNL make interesting to enlarge our 
theoretical knowledge about this last process. 

We have reviewed and studied the vector resonance dominated 0(p 6 ) local amplitudes 
which main role is to determine the low-z region of the invariant di-photon mass spectrum. 
The conclusion is that this contribution is likely to be negligible. The 0(p 6 ) unitarity 
corrections from K + — > 7r + 7r + 7r~ to K + — > 7r + 77 have been computed and found to 
be relevant increasing the branching ratio (for a fixed value of the c amplitude) around 
30 — 40%. The shape of the z-distribution is shown to be sensitive to the evaluated 
corrections (Fig. 2 and Fig. 3). 

In Fig. 1 we have shown the branching ratio for K + — > 7r + 77 as a function of c. 
The included corrections will allow to get a more accurate determination of c once the 
branching ratio of K + — > 7r + 77 is measured and therefore will provide a new independent 
relation between 0{p A ) weak coupling constants that will improve our predictive power 
in this sector. 
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Figure captions 



Fig. 1 : Br(K + — > 7r + 77 ) as a function of the c amplitude. The dashed line cor- 
responds to 0{p A ) xPT amplitude. The full line corresponds to the amplitude including 
the 0{p & ) unitarity corrections. 

Fig. 2 : Comparison of the normalized z-distribution for K + — > 7r + 77 at 0(p 4 ) (dashed 
line) with our 0(p 6 ) correction (full line) for c = 0. 

Fig. 3 : Normalized di-photon mass spectrum of K + — > 7r + 77 for c = —2.3 (naive 
FM, dashed line) and c = (WDM , full line). 

d Br(K + — ► 7r + 77 ) 
Fig. 4 : ^ -g^ '-^ spectrum for c = as a function of \y\. 
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